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ABSTRACT 

Mrk 231, the nearest ( z = 0.0422) quasar, hosts both a galactic-scale wind and a nuclear-scale iron low-ionization 
broad absorption line (FeLoBAL) outflow. We recently obtained a far-ultraviolet (FUV) spectrum of this object 
covering ~ 1 1 50- 1 470 A with the Cosmic Origins Spectrograph on board the Hubble Space Telescope. This spectrum 
is highly peculiar, highlighted by the presence of faint (<2% of predictions based on Her), broad (> 1 0,000 km s _1 
at the base), and highly blueshifted (centroid at ~ —3500 km s _1 ) Lya emission. The FUV continuum emission is 
slightly declining at shorter wavelengths (consistent with F k oc A. 1-7 ) and does not show the presence of any obvious 
photospheric or wind stellar features. Surprisingly, the FUV spectrum also does not show any unambiguous broad 
absorption features. It thus appears to be dominated by the AGN, rather than hot stars, and virtually unfiltered by 
the dusty FeLoBAL screen. The observed Lya emission is best explained if it is produced in the outflowing BAL 
cloud system, while the Balmer lines arise primarily from the standard broad emission line region seen through the 
dusty (Ay ~ 7 mag) broad absorption line region. Two possible geometric models are discussed in the context of 
these new results. 
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1. INTRODUCTION 

At a distance 13 of only 178 Mpc, where 1" = 0.863 kpc, 
Mrk 231 is the nearest quasar known (Boksenberg et al. 
1977). The analysis by Braito et al. (2004) of the joint XMM- 
Newton-BeppoSAX X-ray spectrum of this object has revealed 
the existence of a highly absorbed (Ah ~ 2 x 10 24 cm -2 ), pow- 
erful (intrinsic 2-10 keV luminosity of ~1 x 10 44 erg s _I ) AGN. 
Mrk 231 is also an ultraluminous infrared galaxy (ULIRG) 
with infrared (8-1000 /xm) luminosity log[LiR/L Q ] = 12.54 
and a bolometric luminosity log| Lno\J L©] ~ 12.60. 14 A de- 
tailed analysis of its Spitzer mid-infrared spectrum indicates 
that ~30% of this bolometric luminosity is produced by a cir- 
cumnuclear starburst with SFR ~ 170 M 0 yr _1 , surrounding 
the dominant optical broad-lined quasar (Veilleux et al. 2009). 
The disturbed morphology of the host suggests that both the 
quasar and starburst are triggered by a recent merger event (e.g., 
Hamilton & Keel 1987; Hutchings & Neff 1987; Surace et al. 
1998; Veilleux et al. 2002, 2006). 


13 Based on a redshift z = 0.0422 and a cosmology with Hq = 73 km s 1 
Mpc , Ujnatter — 0.27 . and O vacuum — 0.73. 

14 This estimate for the bolometric luminosity includes all IR and non-IR 
contributions and is derived by simply assuming Lrol = 1.15 Lir, typical for 
ULIRGs (e.g., Sanders & Mirabel 1996). 


In recent years, Mrk 23 1 has become the archetype of galactic- 
scale quasar-driven winds. These outflow events are purported 
to self-regulate the growth of the black hole (BH) and the 
spheroidal component of the galaxy, and explain the relatively 
tight BH-spheroid-mass relation (e.g., Di Matteo et al. 2005; 
Murray et al. 2005; Veilleux et al. 2005 and references therein). 
Evidence in Mrk 23 1 for a powerful kpc-scale wind with outflow 
velocity of up to ~1000 km s _1 and mass outflow rate of the 
order of ~400-1000 M 0 yr “ 1 has now been seen in the neutral 
gas (Rupke et al. 2005; Rupke & Veilleux 2011, 2012), the 
molecular gas (Feruglio et al. 2010; Fischer et al. 2010; Sturm 
et al. 2011; Aalto et al. 2012; Cicone et al. 2012), and the warm 
ionized gas (Lipari et al. 2009; Rupke & Veilleux 201 1, 2012). 
It is not clear whether this galactic-scale wind is physically 
related to the well-known unresolved broad absorption-line 
(BAL) systems in the core of this object (Boksenberg et al. 
1977; Rudy et al. 1985; Hutchings & Neff 1987; Boroson 
et al. 1991; Kollatschny et al. 1992; Forster et al. 1995; Smith 
et al. 1995; Rupke et al. 2002; Gallagher et al. 2002, 2005). 
The Mrk 231 BAL features are detected in several transitions 
including many low-ionization species (e.g., Nai D 7.7.5890, 
5896, He i* A3889, Can H and K, and Mgn A2800, as well as 
Fen UV1 a 25 99 and Fen UV2 7.2382, hence fitting the rare 
FeLoBAL category). They show strong absorption across the 
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range 3500 to —5500 km s _1 and no significant absorption 
at velocities above —3500 km s -1 (see, e.g., Rupke et al. 2002). 

Mrk 23 1 was the first target observed under our far- 
ultraviolet (FUV) spectroscopic survey of the nearby QUEST 
(guasar-C/LIRG Evolution Study; e.g., Veilleux 2012 and ref- 
erences therein) QSOs with the Hubble Space Telescope ( HST ; 
PID 12569, PI: Veilleux). The unprecedented sensitivity of the 
Cosmic Origins Spectrograph (COS; Green et al. 2012) allowed 
us to obtain the first high signal-to-noise (S/N) ratio spectrum 
of Mrk 231 below ~1500 A. This paper describes the results 
from our analysis of this spectrum. The results from the survey 
will be presented in future papers. The acquisition and reduc- 
tion of the COS data and complementary optical spectrum of 
Mrk 23 1 are discussed in Section 2, followed by a description 
of the results (Section 3) and a discussion of the implications 
(Section 4). Throughout this paper, we use z — 0.0422 as the 
rest-frame velocity of Mrk 23 1 ; this value is considered quite 
secure since it is based on rotation curves derived from the ion- 
ized, neutral, and molecular gas components and the peak of the 
Hi 21 cm absorption feature (see Rupke & Veilleux 2011 for 
references). 

2. OBSERVATIONS 
2.1. HST-COS FUV Spectroscopy 

Mrk 231 was observed for a total 12,539 s spread over five 
orbits using the G130M grating and the 2" 5 x 2" 5 aperture to 
cover ~ 1 150-1470 A at a spectral resolution of ~ 1 5-20 km s _1 . 
All of these observations were done in TIME-TAG mode to 
provide temporal sampling and allow us to exclude poor quality 
data and get improved thermal correction and background 
removal. In the end, we did not have to exclude any data for 
Mrk 231. To reduce the fixed pattern noise and fill up the 
wavelength hole produced by the chip gap without excessive 
overheads, we split each orbit into four separate segments of 
similar durations at two different FP_POS offset positions (either 
segments 1 and 3 or segments 2 and 4, depending on the orbit) 
and two different central wavelengths (1309 or 1327 A). All 
observations were carried out on 201 1 October 15 and processed 
using the standard calcos 15 calibration pipeline. Flat fielding, 
alignment, and co-addition of the individual exposures were 
performed using the IDL routines provided by the COS GTO 
team. 16 

2.2. Ground-based Optical Spectroscopy 

The highest-velocity (~— 8000 km s _1 ) component of the 
BAL systems in Mrk 231 has proven to be time variable 
(Hutchings & Neff 1987; Boroson et al. 1991 ; Kollatschny et al. 
1992; Forster et al. 1995; Rupke et al. 2002, 2005), so we decided 
to get a new optical spectrum of this object to complement 
our 2001, 2004, and 2007 spectra published in Rupke et al. 
(2002), Rupke et al. (2005c), and Rupke & Veilleux (2011), 
respectively, and facilitate the task of comparing these data with 
the COS FUV spectrum. The instrumental setup used for these 
ground-based observations was the same as that used for the 
2004 optical spectra: the Ritchey-Chretien spectrograph on the 
May all 4 m telescope at Kitt Peak was used on 2012 April 27 

15 See Chapter 3 of Part II of the HST Data Handbook for COS 
(http://www.stsci.edu/hst/cos/documents/handbooks/datahandbook) for the 
details of the calcos processing. 

16 Routines are available at 

http://casa.colorado.edu/~danforth/science/cos/costools.html, and described in 
Danforth et al. (2010) and Keeney et al. (2012). 


with grating KPC-18C in the first order and a l" 25 slit to obtain 
4 spectra each of 15 minute duration with a spectral resolution 
of 85 km s _1 . Standard IRAF 17 procedures were used to reduce 
and calibrate this spectrum. 

3. RESULTS 

A near-ultraviolet (NUV) COS image of Mrk 231, obtained 
for accurate target acquisition, is consistent with an unresolved 
point source centered on the optical position of the quasar at 
the resolution of the COS imager (FWHM ~ ((.'05). The FUV 
spectrum of Mrk 23 1 is shown in Figures 1 and 2 and is highly 
peculiar. We have been able to identify unambiguously only one 
emission feature intrinsic to Mrk 23 1 : broad, highly blueshifted 
Lya emission (discussed in detail in the last paragraph of this 
section). This lack of FUV features contrasts starkly with the 
rich optical spectrum of this object (shown in Figure 3) and 
with the rest-frame UV spectra of other FeLoBAL quasars 
(e.g.. Hall et al. 2002). These and other quasar spectra indicate 
that Mrk 231 should have strong broad absorption in Lya 
and possibly also emission and absorption from Pv kill 18, 
1128, Sim 7 1206, Nv77 1238, 1243, Sin 71260 (Sin*71265), 
Oi 71302, Sin 71304 (Sin*71309), Cii 71335 (Cii* 71336), 
and Si iv 771394, 1403 within our spectral coverage. However, 
none of these features is clearly detected. The peculiarity of 
the UV spectrum of Mrk 231 was first noted by Hutchings 
& Neff (1987) and Lipari et al. (1994) based on IUE spectra. 
It was further investigated by Smith et al. (1995) using NUV 
FOS spectropolarimetric data and by Gallagher et al. (2002) 
using a short FOS calibration observation in the archives. In 
particular, Gallagher et al. (2002) noted the likely presence of 
broad C iv 771548, 1550 emission affected by a weak BAL, but 
the S /N ratio of the FOS data was insufficient to provide firm 
limits on the equivalent widths of the emission and absorption 
features below ~1450 A. The new COS data show that the 
spectrum of Mrk 231 at <1450 A is also highly unusual. As 
discussed in Section 4.3, we believe that this is the result of 
a complex geometry with a dusty broad absorption line region 
(BALR) obscuring portions of both the continuum source and 
standard broad emission line region (BELR). 

While the COS spectrum of Mrk 23 1 is dominated by Lya 
emission, it does also present a number of faint but significant 
“bumps” and “wiggles.” We are confident that these features are 
real since they are seen in individual exposures with different 
I P POS and central wavelengths during each of the five orbits, 
but they are not seen in the FUV spectra of our other targets, 
which were observed in the same fashion reaching similar S/N 
ratios. Table 1 lists the properties of these weak features along 
with tentative identifications (IDs). None of these IDs is secure 
except for the narrow absorption feature at 1252 A, which is 
most likely a weak narrow Hi absorber at 3530 km s' 1 
associated with either the nuclear outflow in Mrk 23 1 or 1GM 
along our line of sight. The significant width (~2000 km s _1 ) 
of the absorption feature at ~1392 A suggests an AGN origin. 
If this feature forms in the ~-4500 km s” 1 BAL outflow seen 
in the optical Na i, Ca ii, and He i lines, then the rest wavelength 
of the 1392 A feature should be ~ 1 356 A, for which there is no 
plausible ID. The strongest absorption line we expect to find near 
these wavelengths is Cii 71335 (+Cn* 71336), but this implies 
much lower velocities, only +200 (-70) km s _1 , compared to 
the optical BALs. Nonetheless, this ID for the 1392 A feature is 
consistent with the weaker but similar-looking absorption dip at 

17 http://iraf.noao.edu/ 
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Figure 1. HST-C OS FUV spectrum of Mrk 231, binned by 10 spectral pixels to help show the fainter features. It has not been corrected for the small foreground 
Galactic extinction (Ay ~ 0.03 mag). This spectrum is dominated by broad, highly blueshifted Lya emission. The FUV continuum emission is nearly featureless and 
only slowly declining at shorter wavelengths, consistent with F\ oc A 1 7 (shown as a red line in the figure). It is thus dominated by the AGN and only slightly affected 
by dust reddening (Ay ~ 0.5 mag.) from the neutral and dusty broad absorption line region. 

(A color version of this figure is available in the online journal.) 


Velocity Shift from Systemic (km s -1 ) 

-10000 -5000 0 5000 



Figure 2. Lya emission profile in Mrk 23 1 . The vertical dashed line indicates 
the expected position of rest-frame Lya, while the dotted line indicates a 
conservatively high estimate of the continuum level used to produce Figure 4. 
Note the large blueshift of the Lya emission with respect to systemic (the 
centroid of the profile is at 3500 km s -1 ) and the width of the line 
(>10,000 km s _1 at the base of the profile) which indicates an AGN origin. 
There is no indication that Sin A.A. 1 190, 1193, Ni 11200, and Sim 11207 are 
significantly affecting the profile of Lyo'. The feature at restframe wavelength 
1243.5 A is unidentified, but it could plausibly be highly redshifted Lya emission 
at ~+6900 km s -1 (see Section 3 and Table 1 for more detail). 

~ 1452 A, which might be Si iv X 1394 at a similar near-systemic 
velocity. The apparent absence of absorption in the other half 
of the Si iv doublet, at 1403 A, is consistent with its factor of 
two lower /-value and the lower S/N ratio in this part of the 
measured spectrum. The power-law fit in Figure 1 also suggests 
the presence of shallow absorption features around ~ 1 355 A 



Figure 3. UV-optical spectral energy distribution (SED) of Mrk 231 in log kf\ 
vs. log A.(observed) units. This SED combines the new COS data with archival 
FOS spectra and the 2001 Keck/ESI optical spectrum published in Rupke et al. 
(2002). The old pre-COSTAR FOS NUV spectrum at 2500-3300 A was scaled 
up to match the data at shorter wavelengths (see Section 4.1 for more detail). A 
comparison with our new KPNO spectrum shows no evidence for variations in 
the optical continuum and line emission of Mrk 231. The KPNO spectrum is not 
shown here because it does not cover as broad a wavelength range as the Keck 
data. The blue curves indicate the HST /FOS composite quasar spectrum used 
in the HST COS ETC (shown in gray) reddened by Ay = 1.0, 4.0, and 7.0 mag 
using Ry = Ay / E(B — V) = 3.1 and an extinction curve without the 2175 A 
dust absorption feature (Conroy et al. 2010), normalized to match the observed 
spectrum of Mrk 231 at ~3600 A. Note the dramatic drop in flux in Mrk 231 
from ~4000 A to ~2500 A, requiring dust reddening Ay ~ 7 mag, followed 
by a slowly declining continuum below that wavelength, only slightly affected 
by dust (Ay ~ 0.5 mag). This is consistent with the near-zero polarization at 
~1800 A found by Smith et al. (1995). 

(A color version of this figure is available in the online journal.) 
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Table 1 

Uncertain Features in the COS Spectrum 


7obs 

(A) 

(1) 

3- rest 

(A) 

(2) 

EW 

(A) 

(3) 

Gm 

(10~ 15 erg s _1 crrU 2 ) 

(4) 

Possible IDs 
(5) 

Velocity 
(km s _1 ) 
(6) 

Confidence 

(7) 

1181.3 (0.8) 

1133.5 

-3.3 (0.3) 

3.6 (0.8) 

Ni 7.1134.4 

-240 

3 





Feu 7.1 144.0 

-2750 

3 

1225.4 (0.1) 

1175.8 

+ 0.8 (0.2) 


Cm 7.1175.6 

+50 

3 

1252.0 (0.4) 

1201.3 

—29.27— 46.6 a 

49785 a 

Lya 7.1215.7 

-3550 

0 

1252.1 (0.05) 

1201.4 

+ 0.13 (0.02) 


Lya 7.1215.7 

-3530 

1 

1296.0 (0.3) 

1243.5 

-3.6 (0.8) 

6.4 (1.4) 

Lya 7.1215.7 

+6900 

3 





N v 7.1242.8 

+170 

3 





Sin 71260.4 UV4 

-4020 

3 

1324.0 (0.2) 

1270.4 

-2.6 (0.5) 

2.7 (0.5) 

??? 

??? 

3 

1339.5 (0.5) 

1285.2 

-4.1 (1.2) 

4.7 (0.9) 

Oi 71302.2 + Sin 71304.4 

-4200 

3 

1355.2 (0.1) 

1300.3 

+ 2.7 (0.6) 


Oi 71302.2 + Sin 71304.4 

-740 

3 

1375.0 (0.7) 

1319.3 

-8.2 (2.0) 

13.0(1.7) 

Nin 71317.2 

+480 

3 





C 1 71328.8 

-2140 

3 





Cii 71334.5 UV1 

-3410 

3 

1391.8 (0.1) 

1335.4 

+ 1.8 (0.2) 


Cii 71334.5 UV1 

+200 

2 





Cii* 71335.7 UV1 

-70 

2 

1452.0 (0.4) 

1393.2 

+ 0.7 (0.1) 


Si iv 71393.8 

-130 

2 

1462.0(1.0) 

1402.8 

+ 0.2 (0.2) 


Si iv 71402.8 

-0 

2 


Notes. Column 1 : wavelength of the feature (centroid) in the observer’s frame in A. Column 2: wavelength of the feature (centroid) in the rest frame (z = 
0.0422) in A. Column 3: equivalent width in A measured in the observer’s frame (negative value corresponds to emission). Column 4: flux of the emission 
feature in 1 0 1 s erg s -1 cm -2 . Column 5: possible line identification based on, e.g., Morton (2003), Laor et al. ( 1997), Zheng et al. (1997), vanden Berk et al. 
(2001), and Hall et al. (2002). Column 6: velocity of the feature (centroid) corresponding to the tentative line identification. Column 7: level of confidence of 
the line identification: 0 = certain, 1 = confident, 2 = uncertain, 3 = unlikely. See Section 3 for more detail. 

a The first number corresponds to the core of the Lya emission line, while the second number also includes the contribution from the broad wings. 


(~1300 A in the rest frame) and ~1225 A (~ 1176 A; partly 
affected by the damping wings of Galactic Lya absorption). 
This latter feature may be stellar photospheric Chi XI 176 
near systemic velocity; we return to this identification in 
Section 4.1. 

The three emission features at ~1 181, 1296, and 1375 A 
could not be matched unambiguously to the redshifted atomic 
transitions expected from Mrk 23 1 (or the “template” narrow- 
line Seyfert 1; e.g., Laor et al. 1997). The feature at 1 181 A may 
be N i XI 134 blueshifted by ~— 200 km s~ 1 with respect to the 
systemic velocity of Mrk 231, but this identification is doubtful 
since the expected stronger N i X 1200 transition is not clearly de- 
tected (although blending with the highly blueshifted Lya emis- 
sion could hide some N i XI 200 emission). Another possibility is 
Fe II X1144 blueshifted by —2750 km s -1 . The broad feature 
at 1375 A could plausibly be a blend of Nin X1317 and Ci 
XI 329, but with peculiar shifts of about +500 and -2200 km s' 1 
with respect to systemic, respectively. We consider this identi- 
fication unlikely because (1) these lines are not seen in I Zw 
1 and (2) given the ionization potentials and abundances of Ni 
and C, many other lines should also be present in the COS 
spectrum if those IDs were correct, but they are not detected. 
Highly blueshifted (—3410 km s -1 ) Cii X1334.5 is perhaps a 
more likely possibility. The 1296 A feature could conceivably be 
residual N v X1243 (+ highly redshifted Lya) emission affected 
by strong N v BAL, but the N v + Lya + continuum emission 
over ~ 1270-1290 A would have to conspire with the N v BAL 
to produce a nearly featureless pseudo-continuum and explain 
the absence of the N v XI 238 doublet line, which seems very 
unlikely. Highly blueshifted (—4000 km s -1 ) Sin X1260 may 
also contribute to this feature. Similarly, the weaker feature near 
1339 A may be a blend of Oi X1302 and Si n X 1 304 blueshifted 
by —4500 km s -1 . These large blueshifts would be unusual for 
quasars, but as we discuss below, they would be similar to the 


shift measured in the Lya emission feature. A similar blueshift 
(—3000 km s' 1 ) has been reported in C hi X 1 909, but this feature 
is badly blended with Fe in lines at shorter wavelengths (Smith 
et al. 1995; Gallagher et al. 2002). The Civ emission feature in 
the data of Gallagher et al. (2002) also seems blueshifted, but 
the blueshifted C iv absorption feature complicates the analysis. 
We return to this point in Section 4.3. 

The Lya emission feature of Mrk 231 is very unusual 
(Figure 2). It is centered at ~1252 A or 1201 A in the rest- 
frame, i.e., blueshifted by ~3500 km s 1 with respect to rest 
(assuming no contribution from Sin XXI 190, 1 193, N iX1200, 
and SimX1207 emission). It has a very peculiar profile with 
a narrow peak centered at ~ 1261 A (~— 1500 km s -1 in the 
rest frame) and wings extending over at least 1240-1265 A 
(~6000 km s -1 ) and quite possibly as much as 1225-1280 A 
(~ 13,000 km s -1 )- The broad width clearly points to an AGN 
origin for this feature. The Lya emission is centered on the 
NUV nucleus and spatially unresolved with the rather coarse 
spatial resolution of COS at 1250 A for our settings (~1"3; 
Ghavamian et al. 2010). The integrated Lya flux is ~4.9 x 
10 -14 erg s -1 cuT 2 in the line core and 8.5 x 10~ 14 if one 
includes the broad wings (with corresponding rest equivalent 
widths of 29.2 and 46.6 A, respectively; see Table 1). In 
comparison, the Ha flux measured from our 2012 KPNO 
spectrum is ~ 140 x 10 14 ergs -1 cm 2 (the same within the 
uncertainties as that measured from our older ESI spectrum). 
The Lya /Ha flux ratio of Mrk 231 is therefore only ~0.5% 
of the high-density case B value of about 13. Collisional 
suppression of the level 2 population in hydrogen is known 
to be important in the high-density BELRs of quasars but 
reduces the Lya /Ha ratio by at most a factor of ~2-3 relative 
to the Case B value (e.g., Netzer et al. 1995). The integrated 
Lya/Ha ratio in Mrk 231 is thus <2% the values typically 
measured in BLRs. This already small ratio is further reduced 
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Figure 4. Comparison of the continuum-subtracted Ha emission (red), Lya 
emission (black), and Nai D absorption (blue). The conservatively high 
continuum level shown in Figure 2 was used to derive the Lya emission profile 
shown here. The optical features are derived from our most recent KPNO data 
to reduce profile variability. All three features are on the same velocity and 
absolute flux scales. The H« profile was interpolated at ~ 1000-2000 km s to 
remove the telluric O 2 A.6850 band. Ly a shares a stronger resemblance with Na 1 
D absorption than with Ha emission. The emission feature at ~+6900 km s 1 
could conceivably be faint redshifted Lya emission (see Section 3 and Table 1 
for more detail). 

(A color version of this figure is available in the online journal.) 

by at least an order of magnitude if we consider only the 
Ha emission in the core of the line compared to Lya at the 
same velocity (see Figure 4 and discussion in Section 4.2). 
As we discuss in the next section, plausible explanations are Lya 
destruction by dust grains and the possibility that the observed 
Ha and Lya lines form in different regions. 

4. DISCUSSION 

4.1. The FUV Continuum Emission 

Figure 3 shows the SED of Mrk 231 from ~1000 A to 
~ 10,000 A, where we have combined the new COS spectrum 
with archival FOS spectra and our 2001 Keck optical spec- 
trum. Our COS spectrum matches up precisely with the 1996 
November G160H spectrum from Gallagher et al. (2002) with- 
out rescaling; there is therefore no evidence for FUV variabil- 
ity. The older G190/G270H observations from 1992 November 
from Smith et al. (1995) are systematically lower, perhaps by as 
much as a factor of two, than our spectra and those of Gallagher 
et al.; this indicates real FUV variations and/or lost flux in 
the older pre-COSTAR observations. In Figure 3, we rescaled 
the Smith et al. G270H spectrum to join smoothly with the 
Gallagher et al. data. At optical wavelengths, a comparison of 
the Keck optical spectrum shown in the figure with our new 
KPNO spectrum shows no evidence for variations in the optical 
continuum and line emission of Mrk 231. The KPNO spectrum 
is not shown in Figure 3 because it does not cover as broad a 
wavelength range as the Keck data. 

In the following, we assume, based on the results from 
spectroscopic and spectropolarimetric studies (e.g., Thompson 
et al. 1980; Schmidt & Miller 1985; Goodrich & Miller 1994; 
Smith et al. 1995, 2004; Gallagher et al. 2005), that the bulk of 
the optical continuum emission in Mrk 23 1 is produced by the 
AGN. The situation at shorter wavelengths is less clear. As first 
pointed out by Smith et al. (1995), there is a dramatic change in 
the SED slope around ~2400 A. As shown in Figure 3, the shape 


of the FUV-NUV spectrum, especially below ~ 1 400 A, is not 
consistent with any reddening curve, particularly those having 
the 2175 A dust absorption feature (e.g., Conroy et al. 2010). 
The dramatic drop in flux from ~4000 A to ~2500 A requires 
Ay ~ 7 mag (or E(B — V) = Ay / 3.1 ~ 2.3 mag.; see Figure 3), 
but the continuum level shortward of ~2400 A lies well above 
the predictions. Smith et al. (1995) attribute the bulk of the 
<2400 A continuum emission to unreddened light from ~10 5 
hot, young OB stars surrounding the active nucleus. They argue 
that this component dilutes the polarized light from the nucleus, 
resulting in the rapid decline in the degree of polarization they 
measure from ~3200 A to ~ 1800 A. 

Our COS spectrum provides crucial new constraints on 
this issue. Hot O and early B stars are expected to pro- 
duce several stellar photospheric and wind features within the 
rest-frame 1100-1400 A: Pv XX1118, 1128, Siiv XXI 122, 
1128, Cm AA 1 176, 1247, Nv AA1238, 1243, OvA.1371, Fev 
AA 1360—1380, and SiivXX1394, 1403 (e.g., Leitherer et al. 
2001; Robert et al. 2003; Vazquez et al. 2004). The lack of 
any of these features, especially C in X 1 1 76 (EW < 0.8 A), in- 
dicates that hot OB stars cannot contribute more than ~40% 
of the FUV continuum emission. This number is based on the 
C hi XI 176 equivalent width of 2 A measured by Vazquez et al. 
(2004) in NGC 1705-1. This estimate for the stellar contribution 
is a conservatively high value since we are assuming here that 
the weak inflection at ~ 1225 A is entirely due to C m XI 176, i.e., 
we neglect the possibility of kinematic substructures in the Lya 
emission mimicking this inflection. The apparently large width 
(>1000 km s -1 FWZI) of this feature makes this identification 
as a stellar feature very unlikely. Indeed, our measurements of 
Si iv XX 1394, 1403 suggest smaller contributions from hot stars. 

We conclude that the AGN is the dominant source of the FUV 
continuum emission in our data. Both thermal and non-thermal 
AGN processes may in principle contribute to the observed 
FUV continuum emission. But, as discussed in Smith et al. 
(1995), synchrotron radiation from the AGN can safely be ruled 
out based on the fact that the polarization measurements of 
the strong emission lines are nearly the same as that of the 
nearby continuum. We thus favor the scenario where the FUV 
light is thermal emission from a geometrically thin, optically 
thick accretion disk (e.g., Davis et al. 2007 and references 
therein). The slowly declining FUV continuum (consistent with 
F) oc X 1 7 , within the COS wavelength range) and the near-zero 
polarization at ~1800 A (Smith et al. 1995) implies that this 
emission is affected only slightly by dust (A v ~ 0.5 mag). This 
puts strong constraints on the dust distribution as we discuss in 
Section 4.3. 

4.2. Broad and Blueshifted Lya Emission 

The large blueshift of the broad Lya emission in Mrk 231 
(~— 3500 km s 1 ) is unusual among AGNs. Most data on broad 
emission-line shifts are based on the C iv XX 1548, 1550 doublet 
rather than Lya since Civ is less affected by the intervening 
IGM (Lya forest), is free of nearby line emission (N v blends 
with Lya at low resolution), and is accessible from the ground 
starting at lower redshifts than Lya. The mean shift of Civ 
among radio-quiet quasars like Mrk 231 is —800 km s' 1 
(-300 for radio-loud systems; Richards et al. 2011). A shift 
of —3500 km s~ falls well beyond the distribution of Civ 
blueshifts measured by Richards et al. (their Figure 2, which 
is based on ~30,000 SDSS quasars). We come to a similar 
conclusion when using Lya directly, although here the statistics 
are considerably weaker (Kramer & Haiman 2009, based on 


5 


The Astrophysical Journal, 764:15 (10pp), 2013 February 10 


Veilleux et al. 



Figure 5. Simple geometric spherical model of the central pc-scale regions 
of Mrk 231 (not drawn to scale). The BELR (in blue) is modeled as a thick 
shell while the BALR (in gray) is a thin, differentially expanding dusty shell 
with a small hole along the line of sight to the central FUV-emitting accretion 
disk (<0.01 pc, shown in black). This small hole lets only a small fraction 
(~5%-10%) of the FUV continuum source through the shell. BELR clouds 
along the line of sight and hot gas (shown in pink) within the inner radius 
of the BELR, fully ionized by the X-ray source (shown in red), contribute to 
the X-ray measured Compton-thick absorbing column density. In this picture, 
the observed FUV continuum emission is the ~5%-10% that goes through the 
hole unabsorbed. All of the BELR emission, produced on larger scale, is either 
completely absorbed by the shell (FUV resonance lines; e.g., Lya) or attenuated 
(NUV-optical lines). The weak C iv BAL measured by others suggests that the 
hole in the shell is actually filled with dust-free ionized material of lower density 
(higher ionization). Most of the Lya produced in the shell itself is absorbed on 
the way out partly by direct absorption and partly due to photon scattering on 
the grains. The small fraction of the Lya photons that makes it out of the shell 
comes from the part of the shell facing the observer, resulting in the observed 
large blueshift (see Section 4.3 for more detail). 

(A color version of this figure is available in the online journal.) 


symmetric, centered on the systemic velocity of Mrk 23 1 and 
with broad wings extending to at least ±7500 km s~ 1 from 
systemic and a slight excess emission around —5000 km s' 1 
(see also Punsly & Lipari 2005). 

The Lya emission profile is a better match to the absorption 
profile of Nai D than to the emission profile of Ho!. While 
the Lya emission extends to redder wavelengths than the Nai 
absorption, both of these features are highly blueshifted, with 
blue wings extending to at least —8000 km s _1 . This comparison 
suggests that the small fraction of the Lya emission which is 
able to escape from Mrk 231 does so preferentially along the 
flow direction of the FeLoBAL gas (e.g., radially away from the 
central engine; see Hall et al. 2004 for a similar explanation of 
the peculiar U V Fe n spectrum of SDSS J09 1 1 03 .49+444630.4). 
This can happen naturally if the emission comes from an outflow 
region with a significant radial velocity gradient (e.g., Hamann 
et al. 1993). If the outflowing BAL systems were instead in 
front of the source of the Lya emission, then we would expect 
the Lya emission to be suppressed, rather than preferentially 
transmitted, at the blueshifted velocities corresponding to the 
foreground Lya BAL (e.g., Verhamme et al. 2006; Schaerer 
et al. 2011). This is seen for instance in z ~ 3 Lyman break 
galaxies, where the Lya emission is redshifted by approximately 
twice the outflow velocity measured from interstellar absorption 
lines (e.g., Shapley et al. 2003). Moreover, the BALR would 
have to cover the BELR without covering the FUV continuum 
source — otherwise, one would have to tune the BAL optical 
depths to get transmitted fluxes nearly equal to the continuum 
at adjacent wavelengths — but we already know from the optical 
spectra (e.g., Rupke et al. 2002) that the BAL gas covers at least 
80%-90% of the optical continuum source (it is a lower limit 
because light from an old stellar population may also contribute 
to the residual light in the Nai D trough; Smith et al. 1995). To 
avoid these difficulties, we favor a scenario where most, if not 
all, of the observed Lya emission is produced in the outflowing 
material itself and is not part of the standard BELR; this scenario 
is discussed next. 


only ~80 quasars): the median blueshift of Lya in these low- 
redshift unobscured quasars is of the order of —400 km s -1 and, 
according to their Figure 5, a shift of —3500 km s -1 is seen in 
only two other (unnamed) objects. The large Lya blueshift in 
Mrk 23 1 is even more unusual when we consider only FeLoBAL 
QSOs. In these objects, the broad absorption trough removes the 
blueshifted Lya emission, resulting in an Lya emission profile 
that is shifted to the red (e.g., Reichard et al. 2003; Trump et al. 
2006; Gibson et al. 2009; Richards et al. 2011). The only objects 
with blueshifted emission lines similar to Lya in Mrk 23 1 are 
the “PHL 1811 analog” quasars, where Civ is blueshifted by 
at least —6000 km s -1 and tends to have a lower-than-average 
equivalent width (e.g., Leighly et al. 2007a, 2007b; Wu et al. 
2011,2012). 

To gain insight into the origin of this large Lya blueshift, it 
is instructive to compare the Lya emission profile of Mrk 23 1 
with the features detected at optical wavelengths. The results 
from this comparison are summarized in Figure 4, where the 
unsealed continuum-subtracted emission profiles of Lya and 
Ha and absorption profile of Nai D derived from our recent 
optical spectrum are displayed on the same velocity scale (the 
profiles of the other Balmer emission lines are not shown here 
because they are much less reliable due to lower S/N and strong 
Fen emission affecting H/i and H y). The difference between 
the profiles of Ha and Lya is stunning. Within the uncertainties 
of the continuum placement, the profile of the Ha emission is 


4.3. Physical Model 

The data on Mrk 23 1 are complex with different wavelengths 
and features each providing different clues on the geometry 
of the central engine. Figures 5 and 6 show two simple 
pictures that attempt to explain the main features of these data, 
particularly the peculiar Lya emission, the slowly declining and 
nearly unpolarized FUV continuum, and the absence of broad 
absorption features in the FUV despite the presence of a strong 
optical-NUV FeLoBAL outflow. These simple models are no 
doubt simplifications of the reality and may not explain all of 
the details in the observations. In both pictures, the outflowing 
BALR acts as a partial screen for the continuum emission and 
as the source of the observed Lya emission. We discuss these 
pictures separately, although the possibility that Mrk 231 is a 
hybrid of both cannot be ruled out. 

Figure 5 (hereafter referred to as the “spherical geometry” 
or “spherical scenario”) is an attempt to explain the data using 
the simplest geometry possible: the standard BELR is a thick 
spherical shell surrounded by a thin dusty shell representing 
the outflowing BALR. The source of the FUV continuum is an 
accretion disk on the scale of ~few x 10 15 cm (<0.01 pc) with 
arbitrary orientation. If this geometry is to explain BAL QSOs in 
general, particularly their low frequency of occurrence among 
QSOs, then the BALR shell must subtend only a small fraction of 
4 : r str or be present for only a fraction of the QSO lifetime (e.g., 
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Figure 6. Physically motivated geometric disk model of the central pc-scale 
regions of Mrk 231 (not drawn to scale). A strong clumpy outflow (shown in 
gray, with a typical flow line shown in solid black) is present, lifted off the 
accretion disk (in dark gray), and radiatively accelerated by the UV radiation 
field from the UV-emitting accretion disk (in black). Gas too close to the central 
X-ray source (shown in pink and dubbed the “hitchhiking gas” by Murray et al. 
1995) is too ionized and thus too transparent to UV radiation to take part in this 
wind; it does however contribute significantly to the X-ray absorbing column 
density. The portion of the outflow which is closest to the UV and X-ray source 
is hot and free of dust; it is responsible for the weak Civ BAL feature and 
some of the broad line emission, particularly the blueshifted Lya. The densest 
and most radiatively shielded environment at the base of the flow is neutral 
and dusty; it is responsible for the strong Nai BAL seen in the optical. The 
boundaries between the various regions labeled in this figure are gradual and 
probably much more irregular than shown here. Two lines of sight are shown for 
illustration. The top one goes through the hot dust-free (Ay ~ 0 mag) BALR, 
where all of the observed Lya emission, most of the C iv emission, and some 
Balmer line emission are produced. This line of sight is not aligned with the gas 
flow in the far-side dust-free BALR, so very little Lya emission from this region 
makes it out in our direction. The FUV continuum detected along this line of 
sight is seen through the hot dust-free BAL, producing the faint C iv BAL. The 
lower line of sight goes through the patchy dusty BALR (Ay ~ 7 mag) where 
most of the optical continuum passes through the Nai BALR clouds (hence 
it is reddened and partly polarized). The small fraction of the optical and UV 
continuum that does not pass through the Na i BALR clouds is unpolarized; this 
accounts for only a small fraction of observed optical continuum, but perhaps a 
significant fraction of the observed FUV flux (see Section 4.3 for more detail). 
(A color version of this figure is available in the online journal.) 

merger-driven evolutionary scenario of Sanders et al. 1988; see 
also Veilleux et al. 2009; note that the covering factor may 
be much larger among dust-reddened quasars like Mrk 231 
and those selected via a FIRST-2MASS cross-correlation; e.g., 
Urrutia et al. 2009; Glikman et al. 2012). In Mrk 231, this shell 
has a small dust-free "hole” along our line of sight to the FUV 
source that lets 5%-10% of the FUV continuum through the 
shell without being attenuated, but a negligible fraction of the 
BELR emission goes through this hole since the projected area 
of the BELR is much larger than that of the UV source. In this 
picture, all of the Lya emission detected in Mrk 23 1 is produced 
in the dusty BALR shell. We return to this point below. 

Figure 6 (hereafter called the “disk geometry” or “disk sce- 
nario”) is more physically motivated but also more complicated. 
It is inspired by Gallagher et al. (2002; their Figure 12) and 
models where the BELR is comprised of two connected com- 
ponents, the extended accretion disk atmosphere and a radia- 
tion pressure-driven wind farther out (e.g., Murray et al. 1995; 
Murray & Chiang 1995; Proga et al. 2000, 2008; Elvis 2000; 
Chelouche & Netzer 2003; Proga & Kallman 2004; Proga 2007; 
Richards et al. 2011; Kruczek et al. 2011; Richards 2012). In 


these models, the relative emission line intensities from these 
two components depend on the physical conditions in the out- 
flow and therefore, perhaps, the ionizing spectral energy distri- 
bution (SED), as seen by the atmosphere of the accretion disk. 
The BELR and BALR are comprised of clouds/filaments of var- 
ious sizes, optical depths, dust content, and ionization structures 
(e.g., dusty cold core + hot extended halo); only broad trends in 
some key physical properties are shown in Figure 6. The bound- 
aries between the BELR and BALR are gradual and probably 
much more irregular than shown in this figure. As illustrated in 
Figure 6 and described below, the optical-UV spectra seem to 
favor a line of sight to Mrk 23 1 that skirts the upper edge of the 
dusty wind. 

As discussed in Section 4.2, the peculiar blueshifted Lya 
emission in Mrk 231 appears to be direct evidence for line 
emission associated with some of the nuclear outflowing mate- 
rial. The Lya emission from a large covering factor outflowing 
optically thick (in the Lyman continuum) shell is predicted to 
be an order of magnitude or more larger than the flux we ob- 
serve. Thus, only a few percent of the line photons need to 
escape on the blue side of the line to agree with the observa- 
tions. This suggests two different ways that are conveyed by the 
two scenarios shown in Figures 5 and 6. The first, illustrated in 
Figure 5, is a large covering factor dusty shell containing very 
low ionization gas (low ionization parameter) and undergoing 
large differential expansion. The dust in this shell absorbs most 
of the locally emitted Lya photons. The only line photons that 
reach the observer are the ones escaping from the part of the 
shell moving in this direction. The second scenario, illustrated in 
Figure 6, is a dust-free shell with a small covering factor and/or 
a small Lyman continuum optical depth. In this case, the intrin- 
sically produced Lya line is very weak and most of it escapes 
the gas. 

A key issue in Figures 5 and 6 is the location of the dust. As 
mentioned earlier, the spectropolarimetric study of Smith et al. 
(1995) has reported near-zero polarization at ~1800 A, which 
indicates that the FUV continuum emission is not significantly 
affected by the dust that is causing the optical continuum 
polarization. The dust distribution must therefore have holes 
or be patchy, covering only a fraction of the FUV continuum 
source. Given the small size of this source, i.e., the inner 
accretion disk (nominally <0.0 1 pc across), the dust patches 
and “holes” between the patches must themselves be very small 
and, therefore, they are likely to reside near the UV source rather 
than in the extended host galaxy (this is also the conclusions of 
the spectropolarimetric studies using different arguments). 

At optical wavelengths, we know that a very low ionization 
BAL outflow (with substantial Nai, whose ionization potential 
is only 5.1 eV) covers at least 80%-90% of the nuclear optical 
continuum source along our lines of sight (e.g., Rupke et al. 
2002). A few simple experiments with the photoionization code 
CLOUDY indicate that dust is not only likely present in this 
BALR, it seems to be required to shield Nai from the intense 
UV (>5 eV or <2400 A) radiation field of the AGN that would 
otherwise photoionize it to Na ii. So, in both Figures 5 and 6, the 
dusty BALR may be the patchy dust absorber causing the optical 
continuum polarization (see also Rudy et al. 1985 and Goodrich 
& Miller 1994 for independent supportive arguments). If the 
FUV continuum source is also not completely covered by the 
dusty BALR, then the net observational result could be just what 
we observe in Mrk 23 1 : a red partly polarized optical continuum 
where at least 80%-90% of the flux is transmitted through the 
dusty BALR plus a blue nearly unpolarized FUV continuum 
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spectrum dominated by the small fraction of the FUV source 
not covered substantially by the dusty BALR. This explanation 
for the blue FUV spectrum is strongly supported by the absence 
of neutral and low-ionization BALs in the FUV, because the 
observed FUV flux is not transmitted through this part of the 
BALR. 

Both suggested scenarios are not without difficulties. In the 
spherical geometry of Figure 5, one concern is that a single 
hole that transmits only ~5%-10% of the FUV continuum 
emission may not be able to explain the considerable amount of 
unpolarized residual light in the troughs of the Na i D absorption 
(10%-20% of the continuum at ~6000 A, depending on aperture 
size; e.g., Rupke et al. 2002). The optical continuum source 
in a thin accretion disk is predicted to be about an order 
of magnitude larger than the source of the FUV continuum 
(size ~k 4,/3 ), so a significantly smaller fraction of the optical 
continuum should make it out through the single hole. A patchy 
absorber could explain this result if the gaps or holes between 
patches are typically smaller than any continuum source, but 
this patchy absorber would also let ~5%-10% of the BELR 
emission through the shell without being attenuated, which 
seems inconsistent with the small Lya /Ha emission ratio. 
However, there are a number of issues with these arguments. 
First, as mentioned earlier, light from an old stellar population 
may contribute significantly to the residual flux in the Nai 
BAL troughs. This would alleviate the problem by reducing the 
fraction of the optical continuum emission from the accretion 
disk that gets through the BALR. The fact that the fraction of 
residual light in the Nai troughs is lower when using smaller 
apertures is consistent with this idea. The slightly positive 
slope (~A 1 7 ) of the FUV continuum, attributed to reddening 
with Ay ~ 0.5 mag, predicts a larger contribution of (nearly) 
unpolarized light at optical wavelengths than in the FUV, which 
would help explain the optical data. Finally, recent microlensing 
determinations of continuum source sizes in quasars (e.g., 
Blackburne et al. 2011) suggest that the source of continuum 
emission is larger by nearly an order of magnitude than the 
predictions from the standard thin disk models and nearly 
independent of wavelengths (~k 0 - 2 instead of a 4 ' 3 ). The simple 
single-hole model may therefore be able to explain the low- 
polarization FUV continuum, the faint residual light in the Nai 
troughs, and the small Lya /Ha emission ratio all at once. 

A second, perhaps more serious, concern with the spherical 
scenario, where a single hole transmits nearly all of the observed 
FUV continuum emission, is the apparent lack of significant 
UV continuum variability in Mrk 231 over a span of at least 
15 years (based on the 1996 FOS and 201 1 COS data). The hole 
will change shape and eventually move across our line of sight 
to the small FUV continuum source, resulting in large changes 
in the amount of FUV emission reaching the observer. But the 
timescale for these changes depends on several factors (e.g., 
radius of the BALR shell, size of the hole, and density contrast 
between the shell and the hole), which can only be explored with 
a proper theoretical treatment that is beyond the scope of the 
present paper. A third concern with Figure 5 has to do with the 
need for the BALR to redden the NUV-optical continuum by 
A y ~ 7 mag and to produce detectable BAL features in not only 
Nai and Can H and K but also Hei* A3889 (which arises from 
a metastable excited state populated by recombination from 
Hen; Leighly et al. 2011). This suggests the need for neutral 
and moderately ionized columns log Ah > 10 22 cm~ 2 . Again, 
it is hard to tell without a proper theoretical treatment if enough 
Lya photons are able to escape from this thick shell. 


In the disk geometry of Figure 6, the dusty BALR is in the 
densest and most radiatively shielded environment at the base of 
the flow. This inner BALR extinguishes our view of the BELR 
and, therefore, helps to suppress the broad emission lines in 
the FUV. Dust emerging from the disk might further suppress 
emission from the BELR, especially for a strong resonance line 
like Lya. The weak blueshifted Lya emission profile that we do 
observe could be exempt from this BALR extinction if it forms a 
little farther out in the BALR, above the inner dusty region (see 
Figure 6). In this picture, the observed Lya emission comes from 
the near-side portion of hot dust-free BALR, while the Balmer 
lines arise primarily from the standard BELR seen through the 
dusty BALR, in addition to contributions from the dust-free 
BALR and possibly also the dusty BALR (hence explaining 
the large Ha/Lya emission-line ratio and a more symmetric 
Ha profile). The observed lack of UV continuum variability in 
Mrk 231 is a weaker concern in this picture because part of 
the UV continuum emission reaches the observer through the 
dust-free BALR, while the rest comes through the holes (rather 
than a single hole) of the patchy Na i BALR. A potentially more 
serious problem is the implied orientation of the disk. In the 
idealized disk geometry of Figure 6, the absence of redshifted 
Lya emission in Mrk 23 1 (except perhaps for the faint emission 
feature at ~ 1 296 A, see Section 3, Table 1) implies a rather 
edge-on view of the disk with i > 45° so that the far-side 
portion of the dust-free BALR is largely obscured by the dusty 
BALR on the near-side. On the other hand, the results from mm- 
wave/infrared studies on larger few 100 pc scales indicate that 
Mrk 231 has a disk viewed rather face on (i = 10° -20°; e.g., 
Downes & Solomon 1998; Davies et al. 2004). These results 
can be reconciled if this 100 pc scale disk is tilted significantly 
with respect to the accretion disk. The presence of a warp has 
in fact been inferred in the inner ~200 pc by Davies et al. 
(2004) and results on smaller scales (e.g., Carilli et al. 1998; 
Ulvestad et al. 1999; Klockner et al. 2003; i = 56°) seem to 
suggest a more edge-on orientation that provides support for this 
picture. In reality, the boundaries between the dusty and dust- 
free BALRs are probably much more irregular than shown in 
that figure, so the constraint on the orientation may be relaxed. 
Moreover, the “puffiness” of the dusty BALR, which sets the 
covering factor of the FeLoBAL outflow, is poorly constrained 
but likely high in Mrk 231. The derived fraction of LoBALs 
among quasars depends on the methods of selection, ranging 
from only a few percents among UV- or optically selected 
quasars to more than 50% in the FIRST-2MASS red quasars, 
more akin to Mrk 231 (e.g., Urrutia et al. 2009; Glikman et al. 
2012). Also as mentioned above, the line emission from the 
BALR could be mostly in the forward (radial) direction if there 
are sufficient velocity gradients to create this avenue for escape. 
Since our line of sight is not aligned with the gas flow in the far- 
side dust-free BALR, it may be that very little emission from this 
region makes it out in our direction. Finally, it is also possible 
that the line emission from the far-side BELR and dust-free 
BALR is attenuated not only by the foreground dusty BALR 
but also by other intervening gas, including the hitchhiking gas 
which is nearly Compton thick at least along some sight lines 
(Braito et al. 2004; Morabito et al. 2011). 

The weak but non- zero C iv BAL detected by Gallagher et al. 
(2002) is expected in both scenarios. In the disk geometry, it 
happens naturally along the line of sight that goes through the 
hot dust-free BALR. In the spherical geometry of Figure 5, 
this feature may be explained if the dust-free hole is filled 
with ionized gas of lower density (higher ionization) that 
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has presumably leaked out of the shell. The lack of obvious 
blueshifted Siiv Ik 1394, 1403 BAL in the presence of weak 
C iv BAL is also not surprising since Si iv has a lower ionization 
potential and lower elemental abundance (in the Sun) than C iv, 
and it is not uncommon to see C iv BALs accompanied by weak 
or absent Si iv in HiBALs, i.e., the more highly ionized absorbers 
(e.g., Trump et al. 2006; Gibson et al. 2009; Allen et al. 2011). 
The lack of Nv 111238, 1243 with Civ is perhaps harder to 
explain in this context since the ionization potential of N +4 is 
77 eV and the N v transitions are generally present in HiBALs, 
although noticeably weaker than C iv. As discussed in Section 3 
(Table 1), the emission feature at 1296 A may conceivably be 
residual N v (+ redshifted Lya) emission affected by broad N v 
absorption, but this scenario is a bit contrived. Perhaps a more 
plausible explanation is that the C iv HiBAL detected in the FOS 
spectra obtained in 1996 has since disappeared without affecting 
the Nai BAL. This sort of HiBAL variability would not be 
unusual (e.g., Gibson et al. 2008; the 0.5-2 keV and 2-10 keV 
fluxes of Mrk 231 have increased by more than 50% from 2001 
to 201 1 ; Piconcelli et al. 2013). New COS spectroscopy covering 
Lya and Civ AA 1548, 1550 simultaneously would easily settle 
the issue. 

5. SUMMARY AND FINAL REMARKS 

We have recently obtained a high S/N FUV spectrum of 
Mrk 231 covering ~1 100-1450 A as part of our spectroscopic 
survey of nearby QSOs with HST-C OS. The spectrum presents 
faint, broad, and highly blueshifted Lya emission. The FUV 
continuum emission is slowly declining at shorter wavelengths 
and shows no sign of stellar photospheric or wind features and 
no obvious BAL troughs like those seen in the NUV-optical 
spectrum of this FeLoBAL QSO. These data are best explained 
if the FUV continuum emission is produced by the accretion 
disk of the AGN and the observed Lya emission is produced in 
the outflowing BAL cloud systems. The BAL clouds act as a 
patchy screen for the FUV continuum emission and as a polar 
scattering agent to reproduce the distinct optical polarization 
signatures in this object. Many of these results may be explained 
using a simple spherical geometry or a physically motivated 
but more complex disk geometry (Figures 5 and 6). Both of 
these scenarios are not without problems, however. A proper 
theoretical treatment beyond the scope of the present paper will 
be needed to decide which scenario, if any, is favored over 
the other. 

One may wonder why other FeLoBAL QSOs do not show the 
same complex UV-optical spectrum as Mrk 231. The detection 
of the Na i BAL in Mrk 23 1 is an important clue since it requires 
dust. In the disk scenario, Mrk 231 may be an extreme case of 
the classic BAL outflow/dusty torus model of Konigl & Kartje 
(1994) with larger column density, higher mass-loss rate, and 
hence more shielding, than usual. This extra shielding allows 
dust to survive at smaller radii. This facilitates the appearance 
of dust (and Na I) in the BALR, as well as the formation of 
a dusty toroidal wind farther out. This extreme outflow also 
leads to more broad line emission from a wind component of 
the BELR (e.g., Reichard et al. 2003; Richards et al. 2011). 
The amount of covering by this dusty BALR, which is inferred 
from the depths of the Ca ii and Na i troughs (Rupke et al. 2002; 
Smith et al. 1995), is consistent with the amount of line-of- 
sight covering needed to explain the residual FUV continuum, 
roughly 80%-90%. 

Note finally that in both the spherical and disk geometries 
of Figures 5 and 6, the BALR is located within the nuclear 


<1-10 pc scale region rather than on galactic kpc scales. 
The location of the BALR is an important issue since the 
outflow energetics scale with the distance of the BALR from 
the central engine and thus have a direct bearing on the role of 
quasar feedback in the evolution of galaxies. Recently, Faucher- 
Giguere et al. (2012) modeled FeLoBAL outflows as ISM clouds 
(on ~kpc scales) that are shocked, shredded, and accelerated by 
a hot (normal HiBAL) quasar wind. This is inspired by the 
observational work of, e.g., Moe et al. (2009), Bautista et al. 
(2010), and Dunn et al. (2010) indicating that some of those 
flows, or some (narrow line) components of them, do exist at 
>kpc distances. This model likely does not apply to Mrk 231 
for the following reasons : ( 1 ) the N a i BALs of Mrk 231 are not 
spatially resolved on sub-kpc scale (Rupke et al. 2002; Rupke 
& Veilleux 2011, 2012). (2) The Nai BALs varied, at least in 
the highest velocity components. This puts Mrk 231 in line 
with several other BAL variability studies and, in particular, 
with the work by Hall et al. (2011) on a variable FeLoBAL, 
that favors pc-scale origins for the BALs. (3) The detached 
high-velocity Nai troughs in Mrk 231, which imply a lack of 
nuclear absorbers with velocities between " — 3500 and 0 km 
s -1 , are difficult to explain in the models of Faucher-Giguere 
et al. (2012), where the cold gas is accelerated by the quasar 
blast wave starting at the near-zero velocity of the ambient 
ISM. (4) Several spectropolarimetric studies of Mrk 23 1 have 
suggested that the FeLoBAL clouds are causing the optical 
continuum polarization, which is due to dust scattering in the 
polar illumination cones of the circumnuclear region, roughly 
along the pc-scale jet axis of Mrk 231. This result would be 
hard to explain if the FeLoBAL clouds were distributed on a 
galactic scale, where the influence of the central engine is much 
more isotropic (Rupke & Veilleux 2011, 2012). While these 
arguments favor a nuclear origin for the FeLoBAL system in 
Mrk 231, none of them rule out the possibility that this type of 
shocked and accelerated ISM model may apply to the spatially 
resolved galactic-scale outflow in this object (Rupke & Veilleux 

2011, 2012; Faucher-Giguere & Quataert 2012). 
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